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Summary 


A theoretical investigation of the electronic structure of the ozone molecule has been made 
by the molecular orbital method, based on self-consistent linear combinations of atomic orbitals. 
The dissociation energy, the dipole moment of the ground state and the transitions to the lowest 
excited levels have been computed. All computations have been carried out in different approxima- 
tions, in the simplest case considering only 2pz-electrons, in the next approximation considering 
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all the 2p-electrons and in the most accurate treatment considering all 2s- and 2p-electrons. 
It is found that the 2pz- and the 2s,2p-approximations give more reliable results than the 2p- 
approximation. Energy levels, ionization potentials and transition probabilities are rather similar 
according to the 2pz- and the 2s,2p-approximations, but dipole moment and details of the elec- 
tron distribution as given by population analysis are quite different in the different approxima- 
tions. 

The results of the computations in the 2s,2p-approximation agree comparatively well with 
experiments. The excellent agreement between computed and experimental dissociation energy 
is, however, believed to be fictitious, as inclusion of 1s-electrons will possibly change the energy 
values by several electron volts. The computed dipole moment has the expected sign and approxi- 
mately the correct magnitude. Computed electron transitions have been tentatively assigned to 
observed absorption bands. 

Preliminary reports of some results of the present investigation have been published in February 
1955 [18] and communicated at the Molecular Quantum Mechanics Conference in Texas, Decem- 
ber 1955. 


A. Introduction 


The electronic structures and spectra of triatomic molecules have been discussed 
previously by Mulliken [1]. Recently, Walsh [2] has also studied this subject as part — 
of an investigation of polyatomic molecules. Finally, Moffitt [3] has made a detailed 
discussion of the ground state of sulphur dioxide as part of an investigation of the 
sulphur—oxygen bond. All these authors have used general arguments to make it 
plausible that the molecules have the suggested electronic structures. The structures 
proposed by the different authors have many similarities, but they are somewhat 
different with regard to the characteristics of the one-electron orbitals and the general 
scheme of the orbital energies. It seemed therefore to be of some interest to make 
a more detailed computation of the electronic structure of these molecules. The 
aim of the present investigation is to make a purely theoretical determination of 
the electronic structure of the ground state of ozone using no experimental parameters 
except the interatomic distances. A preliminary discussion of the lowest excited levels 
will also be given. 

Recently [4, 5], it has been universally accepted that the three oxygen atoms of 
the ozone molecule form an obtuse, isosceles triangle, the shortest O—O distance 
being 1.278 + 0.003 A and the apex angle 116°49’ + 30’. Thus, the steric structure 
of ozone is very similar to those of NO,, SO,, ClO,, and other dioxides of elements 
from the sixth and seventh groups. As the S and O atoms have similar electronic¢ 
structures, a comparison between SO, and O, suggests itself. The most probable struc- 
tural formulae of these oxides are shown in Fig. 1. According to Moffitt [3] the struc- 
tures I and III should contribute about equally to the structure of SO,. The double 
bonds of IIT originate from a four-valent state of the sulphur atom including 3d 
orbitals. A corresponding valence state of the apex oxygen atom of O, has obviously 
much lower probability. It has therefore been suggested [6] that the structures I and _ 
II (two of each) should give an adequate description of O3. The experimental values 
of bond energies and bond lengths indicate that the S-O bond should be classified 
as: essentially a double bond, whereas the O-O bond in QO, is intermediate between a _ 
single and a double bond. The electric dipole moment of O3, 0.58 Debye, is smaller 
than that of SO,, 1.60 Debye, although the structure III represents a smaller moment 
than I. The experimental values are, however, easily understood, when it is remem- 
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Fig. 1. Structure formulae of molecules XOQ,. 


bered that O is much more electronegative than 8S. The suggested formulae of SO, and 
O, are, therefore, consistent with what is known from experiments about hou the 
Bonds and the dipole moments of the molecules. 

The present paper presents an investigation of the ground state of ozone according 
to the molecular orbital (MO) method in the linear combination of atomic orbitals 
(LCAO) approximation. Section B1 contains a classification of the possible wave 
functions and the choice of the representations of the atomic orbitals and the sym- 
metry orbitals. Section B2 includes a discussion of the previously proposed electron 
configurations and the electron configurations considered in the present investiga- 
tion. Section B3 deals with the choice of suitable atomic orbitals. Section B4 sum- 
marizes the energy calculation from self-consistent field MO’s in the LCAO approxima- 
tion. Section B5 summarizes the population analysis according to Mulliken. Section 
B6 outlines the calculations concerning the excited states. The Sections C present 
the results of the calculations and a discussion with reference to experiments and 
previous investigations. The Appendix D contains some details of the calculations 
and numerical values of some of the more than one thousand basic integrals. 


B. General methods of calculation 


Bl. Orbital representations 


The ozone molecule belongs to the symmetry group C€,,, i.e., the molecule has a 
two-fold symmetry axis (z-axis) and two planes of symmetry containing this axis, 
namely the plane of the molecule (yz-plane) and a perpendicular plane (#z-plane). The 
notation of the different symmetry species according to Herzberg [8] is reproduced 
in Table 1. (In the notation of refs. [1]-[3] 5, and 6, are interchanged as compared 
with Herzberg.) 

As ozone does not show any paramagnetism originating from the electronic struc- 
ture [5, 9], there are only closed shells in the ground state of the molecule, and the 
total electronic eigenfunction ‘’ must belong to the species A,. As usual, Y’ will 
be expressed as products of one-electron functions belonging 2 any of the species 
1, Mp, 6, or by. As the aim is to find a V’ that in the first approximation can be written 


Table 1. Character table of the symmetry group Cy. 


Con i C, (z) Oy (yz) Oy (@2) 
Ay a ‘py +1 +1 +1 
Ay, dp +] +I] —] —l 
Beith +1 -] +1 =] 
Bay Os +1 -—l —] +1 
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Fig. 2. Notations of atoms and of axes and corresponding AO’s for the ozone molecule. 


as a single, antisymmetrized product, the one-electron functions must be molecular 
orbitals (MO’s). Furthermore, the MO’s should be expressed as linear combinations 
of atomic orbitals (LCAO-approximation). 

In the present case, it is assumed that the 1s-electrons are completely localized and, 
consequently, only the AO’s of the valence shell are used in the LCAO MO’s (2s,2p- 
approximation). Since each oxygen atom has six electrons in the L shell, ‘ should 
be composed of nine, doubly filled MO’s. Calculations have also been carried out 
with the assumption that the 2s-orbitals are non-bonding and can be treated sepa- 
rately (2p-approximation). In this case ‘’ has been composed of three MO’s obtained 
from 2s AO’s and six MO’s obtained from 2p AO’s. Finally, some calculations have 
been made in an approximation corresponding to the usual treatment of molecules 
with conjugated double bonds. This is denoted the 2pz-approximation. In this case, 
all the a, and b, electrons (G-electrons) have been supposed to have the same distribu- 
tion as in the valence states of the atoms (see (10) and (11) in Section B2), and only 
the a, and b, MO’s have been varied. 

The notation used for the 2s and 2p AOQ’s is as follows. The 2p AO’s are chosen 
as three equivalent orbitals with mutually perpendicular axes of quantization and 
are denoted x, y, z. The 2s AO’s are denoted s. The AOQ’s of the apex atom, O,, are 
distinguished by the suffix a, and the AO’s of the two equivalent atoms, O, and O,, 
by the suffixes 1 and 2. 

The atomic orbitals can be combined to symmetry orbitals that are eigenfunctions 
of the symmetry operators of C,,. Omitting the normalizing factors in the notation, 
the symmetry orbitals can be written: 


M1: 81 +8, 84, 2, +29, 2a, Yi — Yo} 
D:: 8, — 89, 21 — 2a; Yrt Yo, Yas | l 
by ty +2, Xa} | " 
Bg: By — No. 
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Table 2. Transformation matrix for the two sets of AO’s. The values in parenthesis 
apply to the 2p-approximation. 


$y ay Yn vy Sa Za Ya va 
SAE Alar 0 0 0 
O71 OR COS ase sinc 0 0 
Ty 0 -— sinaw« +cose 0 
Ty 0 ) 0 +1 
Ta -1/V3(+1) -V2/V3(0) | ec mae 
Cal ; -1/V3(0) +1/V6(+1/V2) +1/V2(+1/V2) 0 
Ca2 ~1/V3(0) — +1/Ve(+1/V2) -1/V2(-1/V2) 0 


Ita i) 0 0 tr 


As the MO’s also shall be eigenfunctions of €,,, each MO should be composed only 
of the symmetry orbitals belonging to a single species. The number of MO’s that can 
be obtained from the L shell AO’s is the same as the number of symmetry orbitals, 
namely 5 of the species a,, 4 of b,, 2 of b, and 1 of a,. In the following the MO’s 
will be denoted la,, 2a,, ..., 5a,, 1b,, .... 

For the discussion of bond properties it is more convenient to transform the orbitals 
to another set of basic AO’s. Let o,;, og. and 2, be three equivalent orbitals obtained 
by combination of s,, y, and z, and directed as indicated in Fig. 2. Furthermore, 
let o,, 7, and %, denote 2p orbitals, 0, with its axis in the bond direction, 7, with its 
axis perpendicular to this direction but in the molecular plane, and 7, with its axis 
perpendicular to this plane (cf. Fig. 2). In the 2p-approximation 7, is replaced by 
8,, and o,, and o,, are equivalent combinations of y, and z,. The matrix of the trans- 
formation from the 2, y, z ... set of AO’s to the o, 7, Z... set is given in Table 2, where 
the matrix elements in parenthesis apply to the 2p-approximation. 

The symmetry orbitals of the second representation are: 


Ay: 8, +8, 61+ Fg, Fait Caz, 1+, = 


Opes a Seo Our Og, Dadi Caen Tas: 


— 
bo 
~— 


by: iy + 7U9; Ia; | 


My: 71 — 72. 


Sometimes it is desirable to distinguish between different molecular orbitals 
belonging to the same species. Then the symbol of the correlated orbital in a linear 
molecule can be added [1], [2], e.g. a, 05. 


B2. Electron configurations 


The following electron configurations, where the orbitals are given in order of 
increasing orbital energy, 
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(<1, 0g)? (by 0)? (44 Og)® (By Iu)? (447)? (By Fu)” (g%q)” (4 Fg)” (017%) (3) 


and 


(44 6)? (by Oy)? (4. 0g)® (Bay)? (47)? (By Ou)? (42% p)? (01 7%q)” (44 Mu)? (4) 


have been proposed, the first by Mulliken [1] and the other by Walsh [2], for triatomic 
molecules with an apex angle of 120° and 18 valence electrons. The difference between 
(3) and (4) concerns the order of the two outermost orbitals, a, and 6, and the correla- 
tion of a, to an orbital of a linear molecule. Moreover, the orbital energies estimated 
by the two authors are somewhat different. Omitting the o,, 7, ... symbols and leaving 
open the question of the order of the outermost orbitals, (3) and (4) can be written 


(Lay)? (1b,)? (2a)? (1g)? (aq)? (20)? (Lag)? (4a,)? (30,)”. (5) 


(5) is also the electron configuration of SO, proposed by Moffitt [3], who does not 
discuss the order of the orbital energies at all. 

Except the nine filled MO’s of (5) there are three more MO’s that can be formed 
from the valence shell AO’s. According to [1] and [2] the order of these orbitals 
should be 2b,, 5a,, 4b,. The energy of the 2b, orbital should be very close to the 
energy of the highest filled orbital. Therefore, it seemed to be of some interest to 
ascertain that the proposed general orbital scheme is applicable to O, and should 
not be slightly changed for this molecule. Consequently, not only (5) but also some 
other configurations have been investigated. As the three highest orbitals of (5), 
la,, 4a, and 3b,, are supposed to be antibonding, the configurations obtained by 
replacing anyone of these orbitals by 2b, can be expected to give comparatively low 
molecular energy values and have therefore also been investigated. The treated 
configurations are thus: 


Wy: (Lay)? (1b,)? (2a,)? (2b,)? (3a,)? (10g)? (Lag)? (4a,)2 (3b,)?, (6) 
+ ee (same) (4a,)? (3b,)? (26,)?, (7) 
¥;: +s (la,)? (4a)? (26,)?, (8) 
bt - (la,)? (3b,)? (26,)*. (9) 


For calculations in the 2p-approximation it has been assumed that the three lowest 
orbitals, 1a, 1b, and 2a,, are pure s-orbitals, and that the other six orbitals are pure 
p-orbitals. 

For calculations in the 2p7-approximation, where only the a, and 6, orbitals are 
analyzed, some assumption must be made about the valence states of the oxygen 
atoms. The valence states of O, and O, must be reflections, but the state of O, is 


supposed to have another orientation in space than the other two. The assumed 
electron configurations are: 


0: (24,)? (2p7,)2p0,2p7%, (10). 


Og: (28a)? (2p Ha)® 2p 001 2p Caz. (11) — 
The o-orbitals of (11) are of course pure p-orbitals. 
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B3. Atomic orbitals 


When a calculation including superposition of many configurations is carried out, 
the question of the basic atomic orbitals is of minor importance. However, when the 
number of configurations is limited, as in the present case, the choice of atomic 
orbitals is a very important one. It has been pointed out several times (e.g. by the 
author [10] and by Sponer and Léwdin [11]) that the simple analytical atomic orbi- 
tals of the kind introduced by Slater do not agree very well with the more accurate 
numerical orbitals obtained by the self-consistent field method according to Hartree 
and Fock (SCF-orbitals). Undoubtedly, the SCF-orbitals give better results in cal- 
culations of atomic energy levels than do the analytical orbitals. Whether the same 
is true for molecular calculations is, however, not immediately obvious. It has been 
shown by Sponer and Léwdin [11] that the analytical orbitals and the SCF-orbitals 
will yield about the same final result in the 2p%-approximation for ethylene, although 
the individual integrals are markedly different. However, this result does not exhaust 
the comparison of the two kinds of atomic orbitals. It is a common approximation 
to use the same exponential in the analytical expressions of the Slater 2s and 2p or- 
bitals. Every SCF-computation shows, however, that the 2p orbitals are more ex- 
tended than the 2s orbitals. Consequently, the computation of such quantities as 
dipole moment and degree of hybridization should give different results using the 
simplest Slater orbitals and SCF-orbitals. Furthermore, the greatest error in using 
analytical orbitals does not originate from the electron interaction integrals, mainly 
depending on the extension of the orbitals. In fact, the integrals representing the 
nuclear-electronic interaction and the kinetic energy are more critical, since they 
depend mainly on the form of the orbitals in the immediate vicinity of the nuclei. 
These later integrals cancel in the investigation [11], which was not concerned with 
the energy of the ground state. It is true that some energy terms of this kind will be 
subtracted from the molecular energy expression when the binding energy is com- 
puted, but these terms are significant, e.g. for the determination of the degree of 
s,p-hybridization. 

In the present investigation the best numerical SCF-orbitals of the oxygen atom 
[12] have been used. The results ommunicated here will in a future publication be 
compared with results obtained by use of simple, analytical AO’s. 


B4. Energy calculations 


The molecular orbitals have been computed by a self-consistent field procedure. 
This kind of molecular orbital calculation has recently been reviewed, improved and 
partly reformulated by Roothaan [13]. The calculation procedure of the present 
paper follows in principal features that of Roothaan. 

Neglecting the kinetic energy of the nuclei, the Hamiltonian operator of the mole- 
cule, expressed in atomic units, is: 


t t,u »m 
where the one-electron operator U, is the sum of the kinetic energy of the electron 
and of its potential energy in the field of the nuclei and the inner shell electrons. 


The second sum of (12) is the inter-electronic repulsion operator, Vy, = ris Tu being 
the distance between electron ¢ and electron wu. The last sum of (12) is the inter-nuclear 
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repulsion operator. This term can be dropped in the intermediate calculations and 
reconsidered only in the final computation of the molecular energy. 

In the Hartree-Fock approximation, where superposition of configurations is 
neglected, the eigenfunction ‘ of H can be written as one Slater determinant of 
MO’s 9;(t) =|¢>. In the LCAO approximation each MO can be written as a combina- 
tion of symmetry orbitals ¢,(t) =|4>, which are themselves combinations of AO’s 
x, (t) =|r>, (ef. Section B1). Thus 


pi= > bu <ul dd =D ae <li>, | 


: ay (13) 
or d= 5 | uld=Z1y els] 


The set of numbers <|?> is the representative of the MO q; in the symmetry orbital 
representation, and the set <7|7> is the representative of g,; in the atomic orbital 
representation. Using a matrix notation the numbers </i> (<r|7>) are the elements 
of a column matrix. In order to obtain the best LCAO MO’s the numbers <|7> 
(<r|i>) should be determined by a variational treatment, as described below. 

The Hartree-Fock equations for the closed-shell ground state are: 


ID {ul Fly |) — a: (ul S| y> [DF =0 (14) 


with similar expressions in the AO representation. S is the quadratic overlap matrix 
with the elements 


<u|S|v>=<ulyy =f cide. 


F is a quadratic matrix representing the one-electron Hartree-Fock operator, and 
the parameters ¢; are the eigenvalues of F. The matrix elements of F are 


<u| Flv =<p|U +25 —K|py. (15) 


The operator U of (15) is the same as U;, of (12). The matrix 2 J — K representing the 
electron interaction operator has the following elements: 


<u |2T—K| ry = S42 foe) ot (a) Vin Ss (t) iu) ated ty — 


— | ok (t) pF (a) Vin gi (t) 5, (w) dtd ty} 
a 2 ld|+ {2 (uv |Ax)—(ux|Av)} 


t 


=2,<r|dla> {2(ur|sr)—(ur|sy)}, (16) 
where 
Gud a>= 5 Geli <i], a7) 
<r|d|s)= > <r|t> Gils) (<i|s> =<8|i>*), (18) 
and 
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(abcd) =f wa (t) pe (u) Ven po (t) pa (u) dt: dtu. (19) 


y in (19) stands for any of the orbitals gy, ¢ or y. The sums over i in (16), (17), (18) 
and below are taken over all the filled molecular space-orbitals, and the sums over 
x, 4 and r, s are taken over all pairs of basic orbitals. 

The numbers <r|i> are determined by the variational principle leading to the 
secular determinant equation 


Det |<u| F|»> — e<u|S|»>|=0. (20) 


This equation is of the 12th degree in ¢ in the 2s,2p-approximation and of the 9th 
degree in the 2p-approximation. According to the symmetry (cf. Section B1) (20) 
factorizes into four determinants corresponding to the symmetry species a, 5,, by 
and a,. These determinants are in the 2s,2p-approximation of the 5th, 4th, 2nd and 
Ist order and in the 2p-approximation of the 3rd, 3rd, 2nd and Ist order. As the 
matrix elements <u|F|y> are functions of the <r|i>’s, cf. (16)—(18), an iterative 
procedure must be used for the determination of the <r|7>’s. 

In the 2p-approximation a satisfactory consistency was obtained after 6 or 7 
eycles. In the 2s,2p-approximation the convergence was considerably slower. As 
each cycle was rather time-consuming in this approximation, a comparatively low 
degree of accuracy was accepted as sufficient. In the most accurate computations, 
which were carried out with the wave function (6), the uncertainty in the 2p-approxi- 
mation of the <r|i>’s was at worst 0.003 and of <r|d|s> 0.001. In the 2s,2p-approxi- 
mation the uncertainty of <7|i> was always below 0.03 and with a few exceptions 
below 0.01. The uncertainty of <r|d/s> was in most cases below 0.005 but in a few 
cases as high as 0.01. 

Consider any operator M, which can be written as a sum of one-electron operators, 
M=XM,. The expectation value <M)» of M is a linear function of the numbers 


<r|d|s>: 


<M>={¥* M¥dr=2 > <i|8> <s| Mi|r> <r] =2 D <r|d|s> <3] Me|r>. (21) 


i 8 


One operator of this kind is the electric dipole moment « of the ground state: 


= 2a 2 Deri, (22) 


where all the coordinates z should be measured from a common origin. The second 
sum in (22) is taken over the nuclei. 
Another operator of this kind is the total molecular energy in the Hartree-Fock 


approximation: 
H= Ui +3> (25i— Ki) +3 & Zi Sm Rem, (23) 
t t »m 
<> = Emo 


(for the definition of the one-electron operators J, and K;,, introduced by Dirac, cf. 
(16) above). Alternative expressions of Hy; can be given by introduction of the 
eigenvalues ¢, of F, cf. (14). We have 
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Se=>i|0+2I-K|i= > <r|d|s> [<s|U|r>+<s|2J—K|r)], (24) 
i i r,s . 


Emo — >, Ze Zm Rim= > <r] dl 8) <8|U |r + 2 


=2>a—-— > <r|d|s><s|2J—K|r). (25) 


When superposition of configurations is included, the total wave function is 


Y= 2 C.F. (26) 


Before superposing the ’,’s, each ‘’, has been specified by a determination of the 
numbers <r|i> of (13) as described above. Afterwards the C,’s have been determined 
by the variational principle with fixed <r|i> values. This procedure leads to a 
secular determinant equation of the type 


Detsl deel Haven Hovaks0. (27) 
When all the MO’s of the superposed configurations are orthogonal, the MO’s of 


’,, must be the same as those of ‘, with the exception of at most two MO’s. Thus, 
when x= A: 


J Wiwidr=0, (29) 
[ VIM Vidr+0, (30) 
[PE Vu Vidr +0. (31) 


The present investigation includes only superposition of configurations in the 
2pz-approximation. The configurations used for the ground state are (6) and (7). In 
this case only the integral (31) + 0. The lowest root of equation (27) gives the energy 
of the ground state of the molecule. The next root is likely to approximate the first 
doubly excited energy level of the symmetry 14,. 


B5. Population analysis 


The numbers <r|d|s> of (18) can be considered as elements of a matrix d, which 
is closely connected with the charge and bond order matrix of Coulson, Léwdin 
and others and with Mulliken’s population numbers [7]. In the following, Mulliken’s 
population notation will be used, but the populations are in most cases expressed by 
aid of the <r|d|s>’s instead of by the <r|i>’s, used by Mulliken. 

Mulliken’s subtotal, net atomic population n(r) is: 


m(r)=2 3 (rlt> [ry =2 <r] dr), (32) . 


and his sub-total overlap population n(r, t), which is the overlap population of the 
AO’s x, and x, of two different atoms, is 
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m(r, t)=4 > |r) <r|8|t> <t|i>=4 El d|7) <r| S|, br Hie. (33) 


The gross atomic population N(r) of the AO y, is 
N (r)=2 d Gilr> dr] 8|t <t]= 25 <ld|r> r|S|0p, (34) 


where the sum over ¢ is to be taken over all AO’s including x, itself. The <r|S|r>’s 
are equal to 1, as all orbitals y, € and @ are normalized to 1. In (32)-(34) the sum over 
? should be taken over all filled molecular space-orbitals, i.e. every doubly filled MO 
taken once, as in the preceding section. The total gross atomic population on a single 
atom k is obtained by summing over all AO’s Xr, Of that atom: 


N (k)= p2 N (rx). (35) 


It is also of interest to consider the formal charge Q in an AO or on an atom, which is 
the difference between the population number WN referred to the atomic valence state 
and the population number N referred to the molecule: 


Q(r)=N(r)- Nir), — Q(k) = N(k) — N(R). (36) 


N(r) and N (k) are usually integers. In the case of ozone the N (r)’s are only integers 
in the AO representation (2). NV (y,) and WN (z,) of representation (1) are fractions. The 
overlap population (i) of a molecular orbital q(7) is: 


r+t 


n(i)=2 > airy |S Gd. (37) 


n(z) can be split up in the overlap populations (i; k, m) between all pairs of atoms 
k and m: 


k>m k>m 
n(i)= > n(i; k, m)=2 > 3 <t|re> <re| S| tm><tm|t>. (38) 
k,m 2M Th tm ' 


‘Finally, the total overlap population n of the molecule is: 


n= >n(i)= > n(r,t)= > n(k, m), 
i 


r>t k>m 
r,t k,m 


where 
n (k, m)=2 > n (i; k, m)=4 D2 <tm|d| 14> cst 
i Tho tm 


is the total overlap population between the atoms k and m. 


B6. Excited states 


In the present investigation only a preliminary treatment of the excited states 
has been carried out. As is well known, the solution of the Hartree-Fock equations 
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(14) yields both the filled MO’s of the ground state and the unfilled MO’s obtainable 
from 2s and 2p AO’s and orthogonal to the filled MO’s. The unfilled MO’s are 26), 
5a, and 4b,. The electron configurations of the singly excited states have been obtained 
by promotion of one electron from a filled MO to one of the empty MO’s, leaving all 
the other MO’s unchanged. It could not be expected that energy values obtained 
from a description of the excited states by a single configuration of this kind are 
satisfactory. However, a comparison of the transitions calculated in the 2s,2p-, 2p- 
and 2pz-approximations indicates, how the results of the present calculation are 
likely to be changed by a more rigorous treatment, and makes it possible to correlate 
the electronic transitions to the observed bands of the absorption spectrum. 

The energy 'Z,,, of an excited singlet state obtained by promotion of an electron 
from the MO ¢, to the MO q,, can be obtained from the energy E, of the ground state 
as follows (cf. [13]). Using (24) and (16) the vertical ionization energy of an electron 
can be written 


e=d|U+27—K|) 
=GU|D+E LUG ViD(D- Gl <6] PID} 


= UU|D+ 3 {2Ju— Ku}, (40) 


where the sum over 7 as always should be over all MO’s, which are filled in the ground 
state, including 7 =/. (40) is also valid for the eigenvalues ¢ corresponding to empty 
MO’s, in which case 7 of course must be different from J. Therefore, 


TE im = Ey Clee (Em al ie > Kin) + Jim nt Kin 
= By — 6; Til = Sige ob im: (41) 


The probability for a transition between the ground state E, and the excited state 
E;m 18 proportional to the square of the modulus of 


R= [V3 Dn Vinde. (42) 


Because of the orthogonality of ‘’y and ¥’,,, only one term is left in (42), and the» 
terms corresponding to the second sum of (22) vanish, : 


R= gir gndt= 2 <1 a> <q|¥| 8> <s|m). (43) 
According to Mulliken and Rieke [14] the oscillator strength f is 
f = 1.085 x 10% yap R* R = 3.039 x 10-°y R* R, (44) : 
where » is the wave number in cm~' and R is measured in a.u. 
For the discussion of the wave functions (7), (8) and (9) it is also of interest to. 
calculate the energy of the doubly excited states, where both the electrons from DY 
are promoted to @,,. As the ionization energy of the second electron in Mis =e, —dy, 


we have 
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| ea Ey 7 (2e, —Jx) + 2 (Ep -) 2 Jim pie Kin) a J 
me We ee tee ty td a 4 ig, ee Kim (45) 


mm 


Energy values computed from (45) will be compared with values obtained by cal- 
culation of the best LCAO MO’s corresponding to the electron configurations (7)—(9). 


C. Discussion of the results for ozone 


C1. Observables in the ground state 


To settle the question which of the configurations (6)—-(9) gives the best description 
of the ground state of ozone, the best LCAO MO’s of the 2p-approximation were 
determined for each configuration as described in Section B4. The molecular energy 
E was computed from (25). A comparison between the H-values of the four configura- 
tions is given in the first columns of Table 3. It is seen that ‘’, yields the lowest 
energy, but that the other values are quite close to H,. It can therefore be expected 
that a superposition of the configurations should improve the energy value signifi- 
cantly. To investigate the order of magnitude of such a treatment a superposition of 
the configurations (6) and (7) has been studied in the 2pz-approximation: 


Y =0,¥, +0,¥>. 


Y’, and ‘’, could not be included as they have different configurations of the G-elec- 
trons a, and b,. The solution of (27) gave an H-value, which was 4.23 ev lower than 
H,, and the corresponding values of the coefficients were: Cj = 0.70, C2 = 0.30. 
Nevertheless, the main part of the present investigation has been carried out with 
the assumption that ‘f’, alone gives an adequate description of the ground state of 
ozone. 

To verify that ‘’, is superior to the other functions in the 2s,2p-approximation as 
well, a few cycles of the variation procedure were carried out for each of the four 
configurations. Although consistency was not attained, the calculations showed 
that H, is well below the other energy values in this as well as in the other two approxi- 
mations. 


Table 3. Energy En= | WAY ,dt in ev; Y', is obtained from configurations 
(7), (8) or (9). 


Em —E, computed from 


minimizing of En equation (45) 
2p-appr. 2pz-appr. 2s,2p-appr. 2p-appr. 2pm-appr. 
H,—E, 8.29 5.80 11.00 21.45 5.80 
H,—E, 4,28 —_— 8.02 12.86 _ 
E,-£#, 3.52 — 4.10 10.09 —— 


ee ee ee 
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Table 4. Vertical ionization potentials, —e of (14), in ev (2p-approximation), 
ee | eS Se 
MO of vormoved —e calculated by variation of the MO’s in Max. diff. 
between é’s 
electron Y, Y, Y, Y, 
ee 
2b, ae 8.7 11.8 13.3 4.6 | 
36, 14.4 15.3 — 14.5 0.9 
4a, 17.3 16.0 15.5 ak 1.8 
la, 15.3 — 17.2 17.4 2A 
1b, 19.0 19.2 20.3 21.3 rs 
3a, 28.0 26.0 24.8 26.2 3.2 
2b, 31.0 ol ei 30.4 30.0 1.7 


The energies E,, H, and EH, have also been computed from (45), using the best 
LCAO MO’s of ¥;. These values are listed in the last three columns of Table 3.. 
They differ considerably from the more accurate values of the two preceding columns. . 
Of course, it could not be expected that such a computation should give any H-value > 
lower than E,. It is, however, interesting to note that the differences H,, — Hj,, 
computed in the two different ways, are about proportional to each other in the» 
2p-approximation. Furthermore, both ways of computation give HL, > H,> H, in} 
the 2s,2p- as well as in the 2p-approximation. 

It is illuminant to consider some details of the results from the variations of f 
WV, VL, Vs and V,. Table 4 shows the vertical ionization potentials, i.e. minus the 
orbital energies ¢ of (14). The last column shows that the differences between corre- - 
sponding e-values are of the same magnitude for low and high orbitals. The order of 
the e-values is the same in all cases with one exception, i.e. the reversion of the 
order of the orbital energies of la, and 4a, for ‘’, and ‘V’,. . 

It is especially interesting to compare ¥’, and Vy, since they have the same kind 
of a, and 6, orbitals, which are seen to be the G-orbitals. (As the notations o and z 
are used in the sense of Fig. 2 in the present paper, the notations G and 7 are used 
in the same way as the usual o and z symbols of conjugated molecules.) The difference: 
in orbital energies for ‘’, and W’, is never more than 2 ev. Other characteristics of: 
the orbitals are, however, more different. Table 5 reproduces the gross atomic 
populations N (r) of (34) and the overlap populations (i) of (37) of the G- and the 
x-orbitals of the two configurations. It is seen that the change in the distribution of 
the z-electrons from N (z,) = 2.00 in V, to N (%,) = 0.62 in V’; is accompanied by an 
equally large change in the opposite direction of the G-electron distribution from 
XN (G,) = 1.98 in ¥, to UN (G,) = 3.42 in Fj. 

The overlap populations are also quite different. Adopting the assumption tha 
n(t) is a measure of the bond strength [7] one finds from Table 5 that the G-bon 
are considerably stronger in ‘f’, than in V’,, 3b, being less antibonding in VY, and 4a 
being bonding in Y’, but slightly antibonding in ‘f’,. On the other hand, the strength: 
of the z-bonds, represented in Table 5 by Xm(b,) and n(la,), are reversed compar 
to the strengths of the g-bonds, since the sum of the z-orbitals is seen to be antibond-. 
ing in Y’, and bonding in ‘,. If these results are applicable to unsaturated molecule 
in general, they indicate that the z-electron approximation should give satisfacto 
results with reference to energy values including spectral lines (see also Table 12), but 
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Table 5. Gross atomic populations N(r) of (34) and molecular orbital overlap 
populations n(?) of (37) for G- and %-orbitals from V’, and V’, (2p-approximation). 
(aaa Se ee gS ln 


o-orbitals 
i a a 
vy) | ti 13 Y, ¥, 
N (2;) 0.524 1.348 I n (3a) + 0.208 + 0.202 
N (y,) 1.764 1.663 | n (4a,) ~ 0.034 + 0.103 
=x N (c,) 2.288 3.011 Xn (a) + 0.175 + 0.305 
N (2a) 1.702 0.649 n (2b,) + 0.208 + 0.230 
N (Ya) 1.722 1.328 n (3 ,) — 0.320 — 0.138 
= N (Gq) 3.424 1.978 Xn (b;) — 0.112 + 0.092 
7-orbitals 
| 
N&,) 1.693 1.000 | n (1b) + 0.218 + 0.166 
| n (2 by) — — 0.303 
N Ga) 6615 9.000 | Xn (bs) + 0.218 — 0.137 
| n (1a) — 0.036 = 


that the bond strengths and dipole moment calculated from the electron distribution 
of this approximation could give an incorrect picture of the molecule. 

In the following, the results of the calculations by use of ¥’, of (6) are reported 
and discussed. The values of energies and dipole moment computed in different 
approximations are reproduced in Table 6. Here, Hajs; = 3 Hatom — Emo, Where Emo 
is computed from (25) and Latom = 426.69 ev is the energy of the valence shell of the 
oxygen atom, computed from the SCF AO’s of [12]. The experimental value of Eno 
has been obtained from the experimental dissociation energy and the spectroscopical 
value of Fatom = 433.04 ev. It should be noted that the computed value of the dissocia- 
tion energy #4;,; is much larger in the 2p-approximation than in the 2s,2p-approxima- 
tion. This discrepancy demonstrates clearly the importance of the forced hybridiza- 
tion, discussed by Mulliken [15]. Such a hybridization between the 1s-shell and the 
valence orbitals has been computed for LiH by the author [16] and for N, by Scherr 
[17]. The comparatively good agreement between the experimental value and the 
value of H4;,, in the 2s,2p-approximation will possibly disappear in a more complete 
treatment including 1s-orbitals, although in this case the effect of the forced hybrid- 
ization will of course be much smaller than in the case when the interaction of 2s-orbi- 
tals and 2p-orbitals is neglected, as in the 2p-approximation. 

The dipole moments computed in the different approximations are of the same 
magnitude as the experimental value, but the sign is reversed in the 2p-approxima- 
tion compared to the 2s,2p- and the 2pz-approximations. The latter are supposed to 
have the correct sign, as it is in accordance with the sign to be expected from the 
structure formulae I of Fig. 1. The fact that the ~-value of the 2p7-approximation 
is several times larger than the value of the 2s,2p-approximation is seen to be con- 
sistent with the previous discussion of Table 5. 


38 543 


INGA FISCHER-HJALMARS, The electronic structure of the ozone molecule 


Table 6. Observables of the ground state of ozone: Hy, of (25), Eains = 3 Latom @ 
Eno, & of (22) and e of (14). 

a ae b. i Da A eee enc esas je Siphon SE 

Values calculated from ‘’, in 


Experimental 
Observable es Be z 
28,2p-appr. 2p-appr. 2pm-appr. 
ee ee Oe EE eee 
Exot (ev) — 1305.32 — 1285.13 — 1314.01 — 86.99 
Eaiss (ev) + 6.20 + 5.05 + 33.94 + 40.31 
uu (Debye) 0.58(Og 2) 0.25(Og ) 1.06(Og ) 1.26(Og ) 
ae Se ee ee eee 
€ (ev) for 
la, — 47.32 — 46.56 = 
1b, — 40.96 — 39.23 a 
2a, — 32.63 — 32.98 a 
2b, — 25.21 ~ 31.05 — 
3a, — 22.81 ~ 27.99 a 
lb, — 21.26 — 19.02 — 22.24 
3b, — 15.61 — 14.41 “= 
4a, — 14.82 —17.29 a 
la, — 13.66 — 15.26 — 11.65 
per Qby — 2.00 — 0.02 — 2.24 
5a, + 4.31 — 0.94 -- 
4b, ei By ie By + 1.65 — 


The orbital energies ¢ have in most cases about the same value in the different 
approximations, the only exception being the e-values of 2b, and 3a,. The agreement 
between the 2s,2p- and the 2pz-approximations is especially interesting in view of 
the previous discussion. 

It should, however, be noted that the order of the outermost orbital energies is 
different in the 2p- and the 2s,2p-approximations. As no values of experimental ioni- 
zation energies have been published, it is difficult to tell, which order is the correct 
one. Furthermore, the order could be changed once more, when interaction with the 
1s-electrons is considered. In fact, for N, [17] this interaction increased the computed 
separation of the 30, and the 1z,, orbitals by 0.8 ev. The separation of the orbitals 
of ozone is of this magnitude: Ae(4a, — 3b,) =0.8 ev and Ae(la, — 4a,) = 1.1 ev. 
In addition, several integrals are computed by rather approximate methods, cf. 
Section D4. This lack of accuracy may have different effects on the different symmetry 
species. Judging from the computed e-values it is therefore impossible to give pre- 
ference to any of the configurations (3) or (4). 


C2. +The molecular orbitals in the ground state 


In Table 7 and Table 8 the self-consistent coefficients <u|i> of (13) for the LCAO 
MO’s are compiled. The coefficients are given for the two representations discussed in 
Section B1. When the computed orbitals are to be compared with the configurations 
(3) and (4), representation (2) should be used. As could be expected, none of the 
orbitals a, or b, can be uniquely correlated to any orbital of a linear molecule, i.e. 
they are not exclusively o-orbitals or z-orbitals. However, most of them are essen- 
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Table 7. LCAO coefficients <u|i> of (13) for the molecular orbitals a, and 6). 
a 


2s,2p-approximation 2p-approximation 
MO Representation (1) 

8, +8, Sq panes | -yal| ate, Sa | Z,+2Z, Za | u- ve 
la, + 0.319 | + 0.738 | — 0.058 | + 0.144 | — 0.131 ||+ 0.385 | + 0.790 — — = 
2a, + 0.971 | — 0.685 | + 0.036 | + 0.118 | + 0.018 || + 1.008 | — 0.735 
3a, — 0.382 | — 0.096 | + 0.196 | + 0.802 | — 0.359 — — +0.264 |— 0.451 ]+ 0.761 
4a, + 0.078 | + 0.166 | + 0.958 | — 0.467 | — 0.232 — — + 0.529 |+ 0.777 |+ 0.421 
5a, — 0.313 | + 0.733 | + 0.362 | + 0.551 | + 1.078 || + 0.819 | — 0.534 | — 0.563 

81 — 82 27 2, Ya | YWtY2 $1 — 8, | 21%, | YrtYe 
1b, + 0.738 — 0.178 | + 0.359 | — 0.145 + 1.000 —_— — — 
2b, + 0.687 + 0.563 | — 0.562 | + 0.120 — — 0.263 | + 0.783 | — 0.309 
30, + 0.052 — 0.433 | — 0.033 | + 0.893 — — 0.050 | + 0.486 | + 0.980 
4b, — 0.565 + 0.763 | + 1.037 | + 0.530 — + 1.009 | + 0.545 | + 0.038 

Representation (2) 
Og Ty Og Wy 

8, +8, | 0, +0, |Oa1tCa2| 11+ 2, Ia 8, +8, | 04+ 6, |Oa1t+Ca2|] 1%, +7] (8a) 
la, + 0.319 | + 0.143 | + 0.685 | — 0.011 | + 0.309 || + 0.385 — — — + 0.790 
2a, + 0.971 | — 0.034 | — 0.491 | — 0.024 | — 0.492 ||+ 1.008 — — — — 0.735 
3a, + 0.382 | — 0.214 | — 0.385 |+ 0.342 | + 0.710 — + 0.792 | + 0.451 | — 0.126 — 
4a, — 0.078 | + 0.276 | + 0.134 | + 0.955 | — 0.477 as — 0.629 | + 0.777 | — 0.268 -- 
5a, + 0.313 }+ 1.116 | — 0.916 | — 0.188 | + 0.027 _ — 0.081 | + 0.534 |+ 0.988 — 

Ou Itg Ou Itg 

8; —8, | O1— Og |Oq1~ Fa2 Ty — Mg 8; — 8, | 01-02 |0a1~ Gaz ly — Wp 
1b, + 0.738 | + 0.212 | + 0.359 + 0.082 + 1.000 — — — 
2b, + 0.687 | — 0.371 | — 0.562 — 0.428 — + 0.369 | + 0.783 + 0.074 
30, + 0.052 | — 0.586 | — 0.033 + 0.822 — + 0.843 | — 0.486 — 0.533 
4b, — 0.565 | — 0.827 | + 1.037 — 0.396 — + 0.507 | — 0.545 _ +0.855 


Table 8. LCAO coefficients oh of (13) for the molecular orbitals a, and by. 


MO 2s,2p-approximation 2p-approximation 2pz-approximation 
Ig ee Ig 
ty—@, Le Ups, 

la, + 1.000 | + 1,000 | + 1.000 

Tu Tu It u 

Uy +2 La y+ X, La 

1b, + 0.539 + 0.704 + 0.770 + 0.456 + 0.359 + 0.837 
2b, + 0.893 — 0.768 — 0.702 + 0.937 + 0.978 — 0.621 
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M W 2s, 2p 2p apmt 


; 
Fig. 3. Energies € and other characteristics of the MO’s of ozone according to Mulliken [1], M, 
Walsh [2], W, and different approximations of the present investigation 2s,2p, 2p and 2p7z. 


tially o- or z-orbitals. In Fig. 3 such assignments have been made. In a few cases the 
o-part and the z-part of the orbitals are about equally large. Then both symbols 
have been given with the dominating one in the first place. 

In Fig. 3 are also reproduced the orbital energies and characteristics given by 
Mulliken [1] and Walsh [2]. Mulliken’s values are taken from his general orbital 
scheme. Walsh has not given any energy scale on his scheme. The arbitrary assump- 
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tion has therefore been made that the lowest and the highest filled orbitals of [2] 
have the same energies as those of [1] and the intermediate e-values have been 
determined by assuming a linear energy scale of the figure in [2]. It should be noted 
that according to the authors the orbital schemes of [1] and [2] are not supposed to 
be quantitatively correct. They are, however, reproduced in Fig. 3 in order to show 
the arrangement of the levels. 

Table 7 and Fig. 3 show that the z-parts of the filled a, and b, orbitals are more 
dominant according to the 2s,2p-approximation than according to the 2p-approxima- 
tion. It is supposed that the z-part will be the most important one also in a more 
accurate treatment. Consequently, the 4a, MO should be essentially z, as suggested 
by Walsh [2]. According to the present computations the level 1b, should lie above 
the levels 3a, and 2b,. It is believed that this result is correct, at least when the d- 
orbitals of the apex atom do not contribute essentially to 1b,. According to my opin- 
ion all the s-orbitals of a 90° molecule, la, =s, +5), 2a, =s, and 1b, =8,— 85, 
should have lower energy than any other valence shell orbitals. This point is supported 
by the fact that in the free oxygen atom the energy of the 2s-orbital is 33.7 ev whereas 
that of the 2p-orbital is only 17.2 ev. Moreover, the bonding pé-orbitals 3a, and 26, 
should be lower than the bonding pz-orbital. Although the spacing of the orbitals 
is likely to be different in the 120° molecule, I believe that they should have the 
same order as in the 90° molecule, which is in agreement with both approximations 
of the present computations. 


C3. Population analysis of the ground state 


The overlap populations n(i;k,m) of (38) and n(i) of (37) for the molecular 
orbitals are collected in Table 9. The assumption of a qualitative correlation between 
the overlap population and the bonding power [7] will be adopted in order to discuss 
the results and compare them with the qualitative description of the orbitals given 
by Walsh [2]. According to this author the la, and 10, orbitals should be essentially 
nonbonding lone pairs s, +s, and s, — s,. According to the 2s,2p-approximation in 
the Tables 7 and 9 they could be better described as strong g-bonds formed by 
s,p-hybrids at both O,, O, and O, with la, centered on O, and 16, on the other 
atoms. Because of the large overlap of the 2s-orbitals with the AO’s of neighbouring 
atoms it does not seem to be appropriate to consider the s-electrons as lone pairs, 
although Table 11 shows that the gross atomic populations of s, and s, are close to 
2.0. The bonding properties expected by Walsh for 2a, and 26, have not been found 
for these orbitals but seem to be more characteristic of the lower orbitals. The next 
higher orbitals of the same species, 3a, and 30,, are slightly O,-O, bonding. However, 
3b, is essentially localized to the pair O,, O, and shows many similarities with lag, 
which could be expected from the general orbital schemes of [1] and [2]. For symmetry 
reasons the a, and b, orbitals cannot be complicated by s,p-hybridization as the 
a, and 6b, orbitals and need no comment. 

The values of the 2p-approximation differ considerably from those of the 2s,2p- 
approximation. However, in both approximations the bonding properties are given 
by the following scheme: 


Lit, x by dy by 
0,-0, bonding antibonding antibonding bonding 
0,-Oa antibonding bonding nonbonding bonding 
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Table 9. Overlap populations n(i; k, m) of (38) and 2 (7%) of (37) for molecular 
orbitals and total overlap population » of (39). 


ee 


2s,2p-approximation 2p-approximation 
MO 
n(i; O,, O2) nt; Oy, Oa) n(2) n(t; O,, O2) n(%; O1, Oa) (2) 
la, + 0.0120 + 0.1530 + 0.3180 + 0.0054 +0.1141 + 0.2337 
2a, + 0.0319 — 0.2118 — 0.3916 + 0.0372 — 0.2786 — 0.5199 
3a, — 0.0014 + 0.0276 + 0.0537 + 0.0611 + 0.0736 + 0.2083 
4a, + 0.0391 — 0.1131 — 0.1871 + 0.0273 — 0.0305 — 0.0337 
(5a,) (+ 0.1780)  (— 0.6425) (—1.1070) (+ 0.0546)  (—0.1364) (—0.2182) 
Ld, + 0.082 — 0.144 — 0.207 + 0.131 — 0.121 — 0.112 
1b, — 0.0445 + 0.1467 + 0.2489 — 0.0395 — — 0.0395 
2b, — 0.0170 — 0.0672 + 0.1514 — 0.0146 +0.1112 + 0.2078 
3b, + 0.0999 +0.0018 — 0.0964 — 0.1223 — 0.0987 — 0.3196 
(46,) (—0.1230) (-0.6023) (— 1.3276) (— 0.0366) (—0.1583) (—0.3531) 
xb, — 0.161 + 0.081 + 0.001 — 0.176 + 0.013 — 0.161 
1b, + 0.0097 +0.1071 + 0.2238 + 0.0198 + 0.0992 + 0.2182 
(2b.) (+ 0.0266)  (—0.1936) (—0.3607) (+ 0.0165) (—0.1858) (-— 0.3551) 
Xb, + 0.010 + 0.107 + 0.224 + 0.020 + 0.099 + 0.218 
la, ~ 0.0358 8 — 0.0358 — 0.0358 = — 0.0358 
Lids — 0.036 — — 0.036 — 0.036 — — 0.036 
“i=Lat+xXb} —0.106 +0.044 n=—0.018 — 0.061 — 0.010 n= -—0.081 
2pz-approximation 
lb, + 0.0043 + 0.0849 +0.1741 
(2b.) (+ 0.0320) (—0.8430) (—0.3110) | 
X be + 0.004 + 0.085 + 0.174 | 
la, — 0.0358 — — 0.0358 
Lid, — 0.036 — — 0.036 
Lti=Xat+xUb} —0.032 +0.085 n=+0.138 


This indicates that there should be a b,¢-bond and a 6,%-bond between the apex 
atom and each of the other atoms. According to Table 9 both bonds are weaker 
than ordinary 6- or Z-bonds, where n(O,, O,) is computed to be + 0.2. It is remarkable 
that the total overlap population n is slightly negative except in the 2pz-approxima-. 
tion, and that the total overlap of the s-orbitals in the 2p-approximation is —0.33. 
It seems as if the antibonding orbitals were somewhat favoured, when the bonding 
properties are estimated by population analysis. It is, however, believed that the 
qualitative picture of the bonding properties, given by Table 9, is correct. 

Table 10 shows the overlap populations n (r,t) of (33) between atomic orbitals of 
representation (2). It is seen that in the 2s,2p-approximation there is a strong repul- 
sion between s, and the orbitals of O,, especially o,,. The large overlap integral 
<8,|S|o,;> =0.44 and the almost complete gross atomic population of 8,, N (8) =. 
1.97, give the explanation of this repulsion. Moreover, the strongest a, and 6, bonding: 
is between o, and o,;. From the point of view of overlap populations it is therefore 
justified to neglect the s,p-hybridization of the bonding orbitals at the atoms O, 
and O,, since o, is a pure p-orbital (cf. Table 2). 
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Table 10. Overlap populations n(r, t) of (33) between atomic orbitals of represen- 
tation (2) and n(k, m) of (39) between atoms k and m. 


2s,2p-approximation 2p-approximation 


n(r, t) or 
k, , 
EAR, 3) Da, Phra va, + be Da, Xb, Li=Ua,+ =I, 
N (8, 82) + 0.0438 — 0.0402 + 0.0035 + 0.0426 — 0.0395 + 0.0031 
1 (815 59) 0.0061 +0.0101  +0.0039 a = i= 
mn (845 to) —0.0011  -—0.0080  —0.0091 < se = 
n (O45 83) 0.0061  +0.0101  +0.0039 Ea = es 
n (0;, Op) + 0.0124 — 0.0549 — 0.0425 + 0.0883 — 0.0906 — 0.0023 
nN (04, Te) — 0.0059 —0.0110 — 0.0168 - 0.0021 — 0.0151 — 0.0172 
” (21, 83) 0.0011  —0.0080  —0,0091 ay es .~ 
N (71, Ge) — 0.0059 — 0.0110 — 0.0168 — 0.0021 — 0.0151 — 0.0172 
1 (Ty, Te) + 0.0517 — 0.0483 + 0.0034 + 0.0044 — 0.0162 — 0.0118 
n(O,, Oz) + 0.082 — 0.161 + 0.131 tha 
| Ud, Py PeeWee ee Ua, Xi be 4= 2 d,+ 2b, 
(71, Tp) — 0.0358 + 0.0097 — 0.0261 — 0.0358 + 0.0198 — 0.0160 
n (Oj, Oz) 0.036  +0.010 ~0.036  +0.020 
| Day; 6, Dt=Da,+2Xb, PEO R Diop 2 = a Oy 
n (81, Far) —0.1804 —0.0552  —0.2357 - és a 
n (815 Fae) —0.0035 +0.0011 0.0025 = = == 
n (815 Ta) ~ 0.0017 Pas ~ 0.0017 — 0.1644 = — 0.1644 
N (01, Car) + 0.0918 + 0.1310 + 0.2228 — 0.0307 — 0.0254 — 0.0561 
n (4, Gas) +0.0212  —0.0303  —0.0091 +0.0088  —0.0073  +0.0015 
n (Oy» Ma) — 0.0306 ae — 0.0306 os te ee 
n (7 Fax) 0.0000 —0.0008  —0.0008 +0.0145 0.0182  —0,0037 
n (M4 Gaz) —0.0001 +0.0354 +0.0353 +0.0504 +0.0634 +0.1139 
n (2s Ma) — 0.0408 Bs ~ 0.0408 = as _ 
n (QO 4, Og) — 0.144 + 0.081 SI a | + 0.0138 
Ld, Zi bg. 204 = Lids $21 by Xi a, Xb, Bi=La,t+ Ub, 
1 (Fy Fa) oe +0.1071  +0.1071 oe +0.0992  +0,0992 
n (Oy, Og) is + 0.107 er + 0.099 


In the 2s,2p-approximation a smaller contribution to the bonding is given by the 


interaction of z, and o,. In the 2p-approximation, however, this interaction gives 
the essential part of the bonding, since o, and o,; are antibonding, which is consistent 
with the large N (r) values of these orbitals (Table 11). The large o-populations in the 
2p-approximation can be understood from the fact that the neglected repulsions 
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Table 11. Gross atomic populations N(r) of (34), N(k) of (35) and charges Q (r) 
and Q(k) of (36). 


o_o 


O 2s,2p-approximation 2p-approximation 
A 
or N (r) or N (k) Ot oF N (r) or N (k) Q(r) or 
om, gs av | em) [me ze a | 0m) 
Zh eh 1.003 0.967 1.969 + 0.031 1.000 1.000 2.000 0.000 
- zy 0.905 0.689 1.544 + 0.222 0.388 0.136 0.524 + 1.242 
§ V1 0.290 0.868 1.158 +0.077 || 0.761 1.003 1.764 — 0.530 
8 Ly 1.000 0.397 1.397 — 0.397 1.000 0.693 1.693 — 0.693 
8 8a 1.874 — 1.874 + 0.126 2.000 — 2.000 0.000 
E Za 1.730 _ 1.730 — 0.730 1.702 — 1.702 — 0.702 
x Ya = 1.054 1.054 — 0.054 — 1.722 1.722 — 0.722 
‘ie — 1.206 1.206 + 0.794 — 0.615 0.615 + 1.385 
pe | ee 1.003 0.967 1.969 + 0.031 1.000 1.000 2.000 0.000 
= ea 0.214 0.626 0.840 + 0.160 0.998 0.820 1.818 — 0.818 
8 Ty 0.981 0.881 1.862 + 0.138 0.151 0.320 0.470 + 1.530 
3 Ty 1.000 0.397 1.397 — 0.397 1.000 0.693 1.693 — 0.693 
5 Tey 1.870 = 1.870 — 0.537 2.000 — 7 /2:000 0.000 
$ hes 0.867 0.527 1.394 — 0.060 0.851 0.861 1.712 — 0.712 
. Caz 0.867 0.527 1.394 — 0.060 | 0.851 0.861 1.712 — 0.712 
Ila — 1.206 1.206 + 0.794 | — 0.615 0.615 + 1.385 
\| 
S O; 6.068 — 0.068 5.981 + 0.019 
x Oa 5.863 + 0.137 | 6.039 — 0.039 
8 | 
= 2p7-approximation 
zB 
oa Ly 1.000 0.214 1.214 — 0.214 
Zi, — 1.572 1.572 + 0.428 


(01, 8,) and (6,1, $;) are very much stronger than the neglected repulsions (z,, s,) and 
(Gag, 5;), Which is seen from the overlap integrals in Table 14. This explains also 
the result (cf. Fig. 3 and Table 7) that in the 2p-approximation the a, and b, orbitals” 
are essentially o, and o,, orbitals respectively. ; 

Table 10 shows also that the negative overlap populations are large compared to 
the positive ones. So are for instance n(s,, 2.) (=(s},8,)) in the 2p-approximation 
and (8, 0,;) in the 2s,2p-approximation larger than any of the bonding terms. 
As mentioned above, the total overlap population is negative, although the computed 
binding energy is positive. This emphasizes the qualitative character of the results of 
the population analysis. 

Table 11 gives the gross atomic populations N(r) of (34) and N (k) of (35) as well 
as the charges Q(r) and Q(k) of (36) in e units. Some of the numbers have been dis- 
cussed above, but a few more points should be mentioned. 
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Mulliken [7] has pointed out that calculations on both CO and H,O have given al- 
most the same value, 0.15, for Q (so) of the oxygen atom, indicating that this amount 
of promotion should be characteristic of the oxygen atom in its stable covalent com- 
pounds. According to Table 11 Q(s,) =0.03 and Q(s,) =0.13. Both these values are 
lower than Q (so) for CO and H,0. This can be explained by the facts that the valence 
state of O, is not quite the same as the valence states of most O atoms in stable coval- 
ent compounds, and that the set of AO’s used in the present investigation is not the 
same as that used for CO and H,0. Computations for O, carried out by use of such 
analytical AO’s that have been used for CO and H,O show a much larger amount of 
promotion than that of Table 11. Consequently, the importance of the choice of basic 
AO’s must be remembered in the interpretation of the population values. 

With reference to the discussion of Table 5 it should be noted that the N (zx) 
values of Table 11 are quite different in the different approximations, showing the 
importance of the ¢,z-interaction for the details of the electron distribution. 

The Q(r) values in the 2p-approximation, listed in Table 11, are seen to be surpris- 
ingly large. All the NV (r) in Q(r) of (36) are computed from the electron configurations 
(10) and (11) of the free oxygen atoms. For the 2p-approximation this is perhaps 
not justified in view of the fact that ”(o,, o,,) of Table 10 is negative, which repulsion 
is not consistent with the assumtions in (10). If (10) is substituted by the physically 
less probable configuration 


O,: (281)? (2p0,)? 2p7, 2pm, (46) 


giving repulsion between o, and o,,, we obtain the more reasonable values Q(c,) = 
+ 0.182 and Q(z,) = + 0.530. Nevertheless the remaining Q (r) values are still somewhat 
larger than in the 2s,2p-approximation in agreement with the larger value of the 
dipole moment (cf. Table 6). However, the formal charges on the atoms, Q(O,) and 
@(O,), which of course are independent of the choice of atomic electron configuration, 
are smailer in the 2p-approximation than in the other approximations. This shows 
that it is difficult to draw any quantitative conclusions about the dipole moments 
from the values of the formal charges. The values of the gross atomic populations 
thus support the impression given by the overlap populations that the information 
obtained from the population analysis is qualitatively but not quantitatively correct. 

The results discussed in Sections Cl, C2 and C3 show that the neglects, inherent in 
the 2p-approximation, are so serious that only the values of the orbital energies 
(Table 6) can be used with any confidence. When details about orbitals and electron 
distribution are wanted, the 2pz-approximation could be used to obtain some informa- 
tion about the 7-orbitals. To obtain information about both G- and Z-orbitals the much 
more laborious 2s,2p-approximation must be used. 


C4. Excited states and absorption spectrum 


The experimental investigations of the ozone spectrum have been reviewed by 
Mulliken [1] and Walsh [2] and are summarized in the first columns of Table 12. 

As mentioned in Section B6 an estimate of the energy of the lowest excited states 
of ozone has been made from (41). The oscillator strength f has been computed from 
(44), inserting the theoretically computed y-value. The results are collected in the 
other columns of Table 12. 
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Table 12. The absorption spectrum of ozone. f= oscillator strength of (44). 


a 


Values computed from 


Experimental Trankition 
vel 2s8,2p-appr. 2p-appr. 2pm-appr. | Excited from 
state |ground state 
A AEH f AE f AE f AE f 1A, of (6) 
A ev ev ev ev 
>7500 <1.2 V.w. V2 lOss 4.1 10-8 —- — 1B, 4a,—>2b, 
by? a 10-4 8.2 10-3 Sa 1052 —- — 14) 4a, 5a, 
5730 2.16 
3525-— 3.5 — 8.4 0.02 15.2 0.02 —- — 1B; 3a, 2b, 
—2900 —4.3 | 13.5 10-3 OnmeLOne —- — 1B, 1b, 5a, 
2550 4.86 0.09 9.1 0.6 9.3 0.3 8.9 0.8 1B, la,—2b, 
11.7 0.6 5.8 0.2 —- — 1B, 3b, > 5a, 
<2000 >6.2 13.8 0.1 12.9 0.1 13.5 0.1 1A, 1b, 2b, 
17.6 0.07 16.2 0.09 —_—- — 14, 3a, 5a, 


To estimate the error in the energy values obtained from (41) a comparison can 
be made with Table 3, where energy values obtained from the analogous equation 
(45) can be checked against more accurate values. This shows that the energy values 
of Table 12 may be incorrect by more than 10 ev. Moreover, Table 3 shows a large 
discrepancy between the different approximations, which is also found in Table 12. 
The values of AZ must therefore be considered as very uncertain. 

It is, however, interesting to notice that the 7- transitions la, > 2b, and 1b, 2b, 
are almost equal in all approximations. It can, therefore, be expected that no drastic 
change of these transitions will be obtained by a more accurate calculation. In order 
to verify this, a calculation with superposition of different configurations has been 
carried out in the 2p7-approximation. In this calculation the following configurations 
have been considered for the ground state: j 


‘Yo: Cy (10g)? (lag)? + C, (20,)?(1a,)? + Cs (16,)? (2b,)?, (47) 


and for the excited state 1B,: 
We: Cry (1b,)® 1atg 2bp + O's (2b,)® 1aty 1p. (48) 


By use of (47) and (48) E, is lowered by 4.3 ev and E, by 3.6 ev, so the difference 
EF, — E, is only changed by 0.7 ev. 

Other transitions, reproduced in Table 12, where the results of the 2s,2p- and the 
2p-approximations are very different, could easily change much more drastically, 
when superposition of configurations is included. It is, therefore, difficult to correlate 
the absorption bands with the electron transitions by comparison of experimental 
and theoretical A Z-values. The correlation is, however, facilitated by the f-values, 
which are much more similar in the different approximations than the A E-values. 
The following tentative assignations could therefore be made. The observed, very 
weak, infrared band could represent the 4a,-—> 2b, transition and the weak, visible 
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band the 4a,— 5a, transition. Walsh [2] has suggested that the 4a,—> 2b, transition 
should give rise to one of these two bands, but he expected the 4a, > 5a, transition 
to give a band in the ultraviolet. According to the present results all transitions to 
5a, are likely to give bands of longer wave lengths than expected by Walsh. The 
weak absorption in the near ultraviolet could originate from the 3a,— 2b, or the 
1b, 5a, transition but most unprobably from the la,— 2b, transition suggested 
by Walsh. The last-mentioned transition has a very large f-value in all approxima- 
tions. Therefore, it could possibly correspond to the strong absorption at 2550 A, 
in which the 36, 5a, transition also could be involved. The 1b,-> 2b, and 3a,—> 5a, 
transitions are likely to belong to the far ultraviolet. Finally, all transitions to 4b, 
should give rise to absorptions at much shorter wave lengths than any of the transi- 
tions listed in Table 12. 


D. Appendix: Details of the calculations 


D1. Calculation of molecular orbital coefficients 


The tedious numerical computations of the present investigation consist partly of 
_ the determination of the coefficients of the SCF MO’s and partly of the computation 
of the basic integrals. All the numerical work has been carried out on a desk machine. 

The coefficients of the SCF MO’s were determined by standard methods. The 
secular equations (20) of 5th and 4th order were solved by the method of James and 
Coolidge [19]. 

Scherr [17] has pointed out that in the case of N, the SCF MO’s can be found 
from a few iterations by use of an extrapolation procedure. A similar, labour-saving 
extrapolation cannot be carried out in the case of O,. The reason seems to be that a 
simple relation between the successive changes of the coefficients exists only, when 
the incompletely filled MO’s belong to a single species, as in the case of Ny. In the 
case of O,, however, the incompletely filled MO’s belong to three different species, 
which interact in a very intricate way. 


D2. Notations and parameter values of basic integrals 


In the following some details of the calculations of the basic integrals will be 
mentioned. All values in the Appendix will be given in atomic units. In these units 
the O,—O, distance is 2.411 + 0.006 a.u. and the O,—O, distance is 4.12 + 0.02 a.u. 
[4, 5]. Since in [12] the radial functions P(r) of the SCF AO’s are mostly given at 
r-intervals of 0.2 a.u., it is hardly possible to compute the integrals for other inter- 
atomic distances R than R =0.2 N a.u., where N is an integer. In the present in- 
vestigation the parameters have, therefore, been chosen as follows: 


Rig = Raz = 2.40 a.u. 
Ry, = 4.20 a.u. 


These R-values correspond to an apex angle of 122°. 

As in Section B4, the suffixes k, 1 and m will be used to distinguish the atoms 
O,, k=1, 2, a, and the suffixes ¢ and uw to distinguish the different electrons. r,,, 
for example, is the distance of electron ¢ from atom O,. 
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The 1s AO’s will be denoted K. For the D shell AO’s the notations of Section Bl 
will be used. 


D3. Classification of integrals 


In this section the different kinds of integrals occurring in the calculations of 
molecular energy (23) and dipole moments (22), (42) are classified. 

The one-electron operator U;, of (23) is the sum of the kinetic energy 7; of the elec- 
tron ¢ and the potential energy of the electron in the field of the nuclei and the inner 
electrons: 


UW, = she sae Ont (49) 


The operator U;,, of (49) is defined by 


+ D{2f 7h (u) Venza (u) dtu y(t)—J x3 (w) Vin yp (w) dtu Za(t)}- (50) 


The penetration integral <r|U;1|s> is, therefore, a sum of a one-electron, potential 
field integral, <r| —Z ric |s>, and two-electron integrals of the general type 


(be|df) ={fxs (t) 7a (u) Vin Yo (t) xr (u) ated Tr. (51) 


In (50) Z, =8 in all approximations. The sum over q consists of a single term in the 
2s,2p-approximation, i.e. g = K,, and of two terms in the 2p-approximation, ie. 
q = K,, and s,. In the 2pz%-approximation the sum is extended over the pa AO’s as 
well. There are three pG AO’s of O, or Og, (2p7;,)?, 2po;,, and two of O,, 29641, 29040, 
ef. (10) and (11). 

The electron interaction operator 2J,— K, of (23) and (16) gives rise to two-elec- 
tron integrals of the general type given in (51). 

The dipole moment operators of (22) and (42) give rise to one-electron integrals 
<b|z,|¢>, <b|a,|c> and <b|y,|c>. For the computation of the total molecular moment 
the apex atom QO, is chosen as the origin of the common coordinate system. The 
coordinates of O, and O, are 2%, y;, %, k=1, 2. 

Finally, the lack of orthogonality between the AO’s gives rise to one-electron, 
two-centre overlap integrals <b |S |c>. 

The necessary integrals are thus partly one-electron and partly two-electron inte- 
grals. Furthermore, they can be classified as one-centre, two-centre and three- 
centre integrals. Only members of the two first-mentioned classes can be computed 
with an accuracy, corresponding to the accuracy of the SCF functions for oxygen, 
given by Hartree, Hartree and Swirles [12]. 


D4. Methods of integral computation 


All one-centre integrals have been computed by numerical integration according 
to Simpson’s rule and using the intervals of the radial distance r given in [12]. 
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All one-electron, two-centre integrals and the two-electron, two-centre Coulomb 
integrals (b,,c;,|d,fm) and hybrid integrals (b,c; \d;, fm) have been computed by expand- 
ing the AO’s in spherical harmonics around one of the centres as described by Lund- 
qvist and Léwdin [20]. The operator rj; was expanded in spherical harmonics as 
usual. After integration over the angular variables most terms in the series disap- 
peared. The final numerical integration over r was carried out by use of intervals of 
0.2 a.u. In the ranges, where the intervals given in [12] are larger, the interpolation 
was carried out by use of a function «re °", where « and f were determined so as 
to fit the given points. 

The two-electron, two-centre exchange integrals (b,¢,,|d;,fm) as well as the one- 
electron, three-centre integrals <b,| —Z,rj:'|c,,> and the two-electron, three-centre 
integrals, (b;¢;\d,f) and (b;,¢,\d,f,), were computed by less accurate methods. In 
all these cases the expansion of the AO’s in spherical harmonics according to [20] 
gives rise to slowly converging, doubly infinite series. These integrals were, therefore, 
computed as follows. 

Each AO was expanded in the eigenfunctions K, s, y, z, x ... of the Hartree-Fock 
operator of an oxygen atom. The series should of course in principle include all 
higher AO’s, 3s, 3p, 3d ... as well. As these higher AO’s are not known for oxygen, 
the series was broken off after the L-shell. For the exchange integral (hy h,|h yh.) of 
two hydrogen atoms it was shown that the contribution from higher shells is of 
lower order of magnitude than the first shell contribution. It is reasonable to assume 
a similar result for the case of oxygen. The break off after the L-shell is, however, 
somewhat unsatisfactory, and the accuracy of these groups of integrals is not as 
high as desirable. As there are about six hundred integrals in these groups, it was 
not possible to carry out any still more time-consuming computation. 

The general formula used for the computation of the three-centre integrals 


(b;.€1| def m) is: 
4 (by ¢;| dx fm) 
=> > <ax| S| o> <re| S| fm> (bi Qu |x Pe) + 
ee 


+> > <be| S| a <re| S| fm> (Qer| de re) + 
17k 


+3 <ae|8 |e) du | S| 7m) Or dum fm) + 
+ >. > <bi | S| a> <d|S|1m> (91 ¢1| ?'m fm). (52) 
q1 ™m 


Expressions similar to (52) have been used for the other groups of integrals. (52) is a 
generalization of the formula given by Mulliken [21]. A similar expression has been 
discussed by Ruedenberg [22]. 


D5. Numerical values of some integrals 


From considerations of space, only integrals over atomic orbitals are listed below, 
and of these the numerous two-electron, three-centre integrals are omitted. The 
integrals over symmetry orbitals can be computed from the tables as indicated in 
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Section Bl. Most of the tables give integrals over AO’s of representation (1). The 
reason for this is that it is inconvenient to express the three-centre integrals in 
the AO’s of the less formal representation (2). The self-consistent field computations 
have, therefore, been carried out in representation (1), and the two-centre integrals 
are accordingly given in this representation. Only the overlap integrals are also 
given in representation (2). 

As it can be easily deduced from symmetry, which integrals are equal, or differ 
only with respect to sign, only one integral of such a group is usually given. 


Table 13. Overlap integrals <b|S|c> of atomic orbitals of representation (1) in a.u. 


Xb Xe 


81 8 
$142 
$1 2 
YWYe2 
Y1 22 
% % 
XX 


+ 0.0380 
+ 0.0828 
0 


‘= 0.1129 


0 
+ 0.0345 
+ 0.0345 


<b, |S | ce> | 


Xo Xe 


81 8q 
$1 Ya 
$1 2 
NYa 
Y12%a 
%4%q 
%Xq 


<b, | S| ca> 


++ 0.2707 
+ 0.3068 
+ 0.1697 
— 0.1477 
— 0.1941 
+ 0.0957 
+ 0.2030 


Table 14. Overlap integrals <b| 8|c> of atomic orbitals of representation (2) in a.u. 


Xd Xe 


<b| Se) 


2s,2p-appr. 


+ 0.0380 
+ 0.0725 
~ 0.0401 
+ 0.0946 
—~ 0.0332 
+ 0.0529 
+ 0.0845 
+ 0.4425 
+ 0.0087 
+ 0.0177 
+ 0.4107 
+ 0.0950 
+ 0.1016 
— 0.0030 
— 0.1420 
+ 0.1450 
+ 0.2030 


2p-appr. 


+ 0.0380 
+ 0.0725 
— 0.0401 
+ 0.0946 
— 0.0332 
+ 0.0529 
+ 0.0345 
+ 0.3369 
— 0.0969 
+ 0.2707 
+ 0.2451 
— 0.0705 


- +0.3506 


— 0.0561 
— 0.1951 
+ 0.0000 
+ 0.2030 


sss 
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Table 15. One-electron potential energy integrals <b|—Z, 177} |c> of atomic orbitals 
of representation (1) in a.u. 


ee eee ee ee) ee 
Xb Xe |-8/nelo> | <6|-8/ratl &> | <6|-8/rae|e> 
ee a eer On) ee ee Pe 

818; — 10.1136 SPS PAM | — 1.9038 
YY — 8.8874 — 3.5503 — 1.9682 

Rg By — 8.8874 — 3.2416 — 1.8625 
U2, — 8.8874 — 3.1054 — 1.8625 
8191 0 + 0.7878 + 0.3032 
832 0 + 0.4359 0 

Yy 2% 0 — 0.2462 0 

8,85 — 0.1466 — 0.1915 — 0.1466 
84 Ya — 0.4365 — 0.3440 — 0.2456 
Buea 0 + 0.0260 0 

Yi Y2 + 0.4472 + 0.4666 + 0.4472 
Y, 2, 0 — 0.0363 0 

Zea — 0.1115 — 0.1119 —0.1115 
Gifs — 0.1115 — 0.1072 =O.1L115 
81 8q — 1.4841 — 1.4841 — 0.6699 

81 Ya — 2.0182 — 1.4376 — 0.7306 
Bien — 1.1166 — 0.7954 — 0.4878 
1 8a + 1.4376 + 2.0182 + 0.8367 
"Ya +° 1.0751 + 1.0751 + 0.3084 

Ua 2a + l.1114 Salli id + 0.5811 

Ba org + 0.7954 + 1.1166 + 0.3426 

2 Yo +t 1.1114 +1.1114 + 0.4693 
Ata — 0.3187 — 0.3187 — 0.2305 
% Xa — 0.9336 — 0.9336 — 0.5043 


Table 16. Kinetic energy integrals <b|7'|c> of atomic orbitals of representation 


(1) in a.u. 
Xo Xo <b| Te) Xo te <b| 2 |e> 
8k Sie + 3.1060 Pk Pk + 2.5422 
81 8 — 0.0163 81 8q + 0.0722 
81 Y2 + 0.0181 81 Ya + 0.1195 
812s 0 84 2a, + 0.0661 
Yi Y2 — 0.0445 YWYVa — 0.1926 
Y1 22 0 Y12a — 0.1459 
2 Ze — 0.0011 224 — 0.0096 
bits — 0.0011 Coy + 0.0712 
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Table 17. One-electron dipole integrals <b|z|c> of atomic orbitals of representa- 
tion (1) in a.u. 


2, =%,= +1.1619 a.u. 


a 


ho te <b| a] ¢> Xo te <b|ale> 
yaw al + 1.1619 Xa Xa 0 

8 24 + 0.6693 8a 2a + 0.6693 
81 8» + 0.0442 81 Sq + 0.1575 
81 Yo + 0.0962 81 Ya + 0.1951 
81 2 + 0.0327 81 2 + 0.2954 
Yy 82 — 0.0962 Y1 Sa — 0.1686 
Yr Yo — 0.1312 Y1Ya — 0.0860 
Ys 2s — 0.0717 Yr 2a — 0.3269 
ase + 0.0327 2 Sg + 0.0964 
21 Yo + 0.0717 21 Ya + 0.1006 
2 29 + 0.0401 2 2a + 0.0555 
ee + 0.0401 XL Ly + 0.1182 


Table 18. One-electron dipole integrals <b| y:|e> of atomic orbitals of representa- 
tion (1) in a.u. 


Yi. = —Y2= +2.1000 a.u. 


Xo Xe | <b| yz |e> 
M11 + 2.1000 
8141 + 0.6693 
8a Ya + 0.6693 
$142 + 0.0669 
$1 8q + 0.2847 
$1 Ya + 0.5401 ; 
$1 2%q + 0.1951 
41 Sa — 0.1150 
WYa — 0.1567 
Y1%a — 0.0868 
2 8 — 0.1686 
“1 Va — 0.3230 
2% + 0.1002 
LX + 0.2136 
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Table 19. One-electron dipole integrals <b|a|c> of atomic orbitals of represen- 
tation (1) in a.u. 


X,=2,=0. 

Xo Xe <b | rg | ¢> 
7 + 0.6693 
81% + 0.0327 
YX — 0.0717 
2 Xo 0 

81 2a, + 0.1874 
Yq — 0.2135 
Pilg — 0.1181 


Table 20. Two-electron, one-centre integrals ( (bi ¢x| di fr) fx) of ane orbitals of 
representation (1) in a.u. 


Xo Xe Xa ws | (by. Ci | Uc fre) 
KEK KK + 4.7473 
KE 88 + 1.1345 
KK yy + 1.0988 
88 8s + 0.7972 
8s YY + 0.7733 
yy yy + 0.8080 
yy zz + 0.7275 
Ks Ks + 0.0777 
Ky Ky + 0.0351 
sy sy + 0.1572 
Ye Yz + 0.0403 
KEK Ks — 0.4798 
Ks 88 — 0.0215 
Ks yy — 0.0151 
Ky sy + 0.0378 
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Table 21. Two-electron, two-centre, Coulomb integrals (bx Cr | Am fm) of atomic or- 
bitals of representation (1) in a.u. 


eS _—— 


Lo Ke La Xf (bs ¢; | de fe) Xo Xe Xa Xs (6, ¢1 | da fa) 


Keke 8985 + 0.2381 K,K, 8a 8a + 0.4159 
K, kK, YoYo + 0.2477 K, Ky Yaa + 0.4435 
K,K, 292 + 0.2329 K, Ky Zata + 0.4052 
Ky VeLq + 0.2329 K, Ky tala + 0.3883 
K,K, $2 Yo + 0.0379 K,K, SaYa + 0.0981 
KK, 8229 0 K,K, Sata + 0.0543 
K,k, Yo 2s 0 K, kK, Yaza + 0.0306 
818, Sy 85 + 0.2412 8181 Sa Sa + 0.4055 
818; YoYo + 0.2471 818; Yaa + 0.4261 
818 ata + 0.2327 81 8; ace + 0.3979 
8181 Ly Wy + 0.2327 8181 Ua + 0.3855 
Br8, $2 Ye + 0.0377 8181 8a Ya + 0.0871 
8184 8329 0 81 8y 8aZza + 0.0482 
8181 Yo2% 0 818, Ya2a + 0.0225 
WV Y2Y2 + 0.2575 Wr YaYa + 0.4379 
Wr aes + 0.2415 Yiy1 ata + 0.4192 
YY Usha + 0.2415 WY Va%a + 0.4004 
YWrvY1 82Yo + 0.0420 WY SaYa + 0.0854 
WY 8229 0 WY Sa%a + 0.0613 
WY" Y2%2 0 WY Ya%a + 0.0265 
2% 2a 2 + 0.2283 2121 2q%a + 0.3840 
224 L_Xe + 0.2276 2124 Lave + 0.3758 
Bie S2Y2 + 0.0352 22 SaYa + 0.0877 
22 8325 0 By 2 82a + 0.0345 
2% Yo%e 0 4% Ya%a + 0.0165 
“ee Loy + 0.2283 UX Siete + 0.3681 
Dia, 82 Yo + 0.0352 LyX, SaYa + 0.0776 
U2, 822 0 Hy Ly 8a2q + 0.0429 
© Ly Y2% 0 UX, Yara + 0.0196 
81Y4 82 Yo — 0.0119 81 Y4 SaYa — 0.0222 
8141 822 0 SY ned — 0.0276 
SY Y22 0 8141 Yara — 0.0113 
8424 8229 + 0.0060 8124 sien + 0.0124 
8421 Yo2%e + 0.0025 8421 Yara + 0.0025 
Ys 24 Yo 2 — 0.0012 Yr 2 Yatra + 0.0007 
81 Xy S_Xy + 0.0060 8,2 ey ta + 0.0277 
80, Yoke + 0.0025 8,2, Yata + 0.0127 
8X ZX, 0 8,2; Zqtq + 0.0070 
YX, YX — 0.0012 YX, Ya%a — 0.0041 
YyXy 2_Ho 0 YX, ZqXq — 0.0058 
2,2, eatee + 0.0007 PES Zao, + 0.0032 
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Table 22. Two-electron, two-centre, exchange integrals (bi Cm | dix fm) of atomic or- 
bitals of representation (1) in a.u. 


Xb Xe tats (by ¢2 | ds fa) Xb Ke a Xs (bsca| dy fa) 
KG 8, nS + 0.0000 Ky 84 Ky 8 + 0.0006 
Ky Y2 Kyi Y2 + 0.0000 Kiya Ky Ya + 0.0012 
Ky 2. Ky 2. + 0.0000 Ky % Gp + 0.0004 
Ky, x, Ky, 2%, + 0.0000 1G oi LGaa + 0.0001 
Ky 8, Kiy2 + 0.0001 TG Kya — 0.0001 
Ky 8, Ky 2 0 Ky sq Ky2q — 0.0001 
Ky Y2 Ki 22 0 Ky Ya Ki%a + 0.0006 
8185 8185 + 0.0014 81 8q 8184 + 0.0654 
$142 81Y2 + 0.0040 $1Ya $1Va + 0.0739 
8129 - 81 2 + 0.0003 812g 832g + 0.0294 
8X 81H. + 0.0003 8, Xq 8, Xq + 0.0098 
81 8 81Y2 + 0.0023 8184 81 Ya + 0.0667 
8189 812 0 8184 81 2q + 0.0369 
81 Yo 812 0 81 Ya 812 + 0.0355 
WY WYa + 0.0078 YiYa YWiYa + 0.0307 
Y1%2 Y1% + 0.0008 Y12a Y1%a + 0.0365 
YX Yi + 0.0008 Yiva Yivq + 0.0097 
Y182 YWr1Y2 + 0.0057 Y18q YiYa + 0.0422 
Yi82 Y122 0 Yi8a Yi%a + 0.0463 
YWY2 Y1%2 0 YiVa WN% + 0.0236 
Pee Ree + 0.0006 eee 2420 + 0.0094 
CEP Sy ha + 0.0000 21 Xq Gin + 0.0036 
% 83 Z1Ys + 0.0003 28a 21 Ya + 0.0303 
248 Pies 0 2 8q 21 2a — 0.0061 
21 Ye 212 0 21Va 21%, — 0.0067 
T,X, ©1X>5 + 0.0006 ©, Xq Xa + 0.0242 
1 Sq U1Yo + 0.0003 0 84 TY + 0.0068 
1 85 X12 0 1 8q Dien + 0.0038 
XyYo X12, 0 21 Ya X14 2q + 0.0048 
8185 YrY2 — 0.0030 818q YWiVa — 0.0406 
S1Y2 Yi 8 — 0.0041 8,Ya Y18q — 0.0681 
8182 Y4 2% 0 81 8q Yiza — 0.0418 
8122 Y1 82 0 8124 Y18q — 0.0411 
S81 Ye Yi22 0 81 Ya YZ — 0.0399 
$122 YrY2 0 812 Y: Ya — 0.0222 
818, 22 + 0.0007 81 8a @1 2a + 0.0118 
812 218 + 0.0001 8124 21 8q — 0.0166 
81Y2 212 + 0.0014 81 Va 2124 + 0.0132 
812 21Ys + 0.0002 8124 21Ya — 0.0276 
YrY2 22 — 0.0019 YiYa 2124 — 0.0085 
Y1 22 Z1Y2 — 0.0001 Y12a 21Ya + 0.0214 
8182 XX + 0.0007 $18, 22a + 0.0350 
8,2, 483 + 0.0001 81% 118 + 0.0061 
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Table 22 (continued) : 
See eee ee ee ee 


Xb Xc ars” (b; Cz | d, fe) | Xb Xe Xa Xf (61 a | d, fa) 
81 Ye XH, + 0.0014 81Ya X1Xq + 0.0367 
8X X1Yo + 0.0002 8 Xq, Ya + 0.0046 
Size Hy Xo 0 8124 Uy + 0.0203 
81 Xo X42 0 81 Xq X42 + 0.0026 
YiY2 Dies — 0.0019 Y1Ya @4Xq — 0.0197 
Yy Lo L1Yo — 0.0001 Yi%a Tia — 0.0019 
Y1 2 Hq 0 Y12a #1 2q — 0.0234 
Yy Xe X12 0 Y1Xq Xy2q — 0.0016 
212 2, Xo + 0.0006 Z42q 24 Xq + 0.0096 
2X 2% + 0.0000 2% 2 2q + 0.0001 


Table 23. Two electron, two-centre hybrid integrals (b;, ¢;|d, fm) of atomic orbitals 
of representation (1) in a.u. 


Xo Xe Xa Xf (61% | d; fa) Xo Xe Xa Xf (b, e|dy fa) 
K,K, 8189 + 0.0187 K,K, 818q + 0.1949 
K, Ky, 81Y2 + 0.0535 K,K, 81 Ya + 0.2657 
K,K, 8125 0 K, kK, 82a + 0.1470 
K,K, Y1 82 — 0.0309 K,K, Yi 8a — 0.1869 
K,K, YY — 0.0551 K,K, Y1¥a — 0.1393 
K,K, Yi 2 0 K,K, 12a — 0.1442 
K,K, 2183 0 K,K, 28a — 0.1034 
K,K, 2142 0 K,K, 21Ya — 0.1442 
K,K, 2 Zs + 0.0140 K,K, 2a » +0.0415 
K,K, © 25 + 0.0140 K,K, U2 + 0.1213 
8184 8182 + 0.0174 818, 8184 + 0.1665 
818, 8142 + 0.0466 818, 81Ya + 0.2189 
813; 8129 0 818; 812q +0.1211 
818, Y18 — 0.0300 8,84 Y18q — 0.1644 
818, YWsY2 — 0.0519 8184 Ya — 0.1150 
$181 Y1%2 0 81 81 Y12a — 0.1228 
818, 2489 0 8,8, 218¢ — 0.0910 
818; Z1Ys 0 887 21Ya — 0.1228 
818, 229 + 0.0133 818, 2124 + 0.0390 
818, %1X_ + 0.0133 818, %2q + 0.1070 
WY 818, +0.0179 YWrY1 81 8q + 0.1664 
WM 8142 + 0.0469 WY 81Ya + 0.2145 
WY 812 0 Wi 81 2q + 0.1244 
WM Yi8s — 0.0314 WV Yi 8a — 0.1717 
Ww WYe — 0.0555 wr "Ya — 0.1185 
WY Yi 22 0 WY Yi2q — 0.1318 
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Table 23 (continued) 
eee eee ee ee eee ee 


Xo Xe Xaxs (b,c, | d; fa) 10 te tats (b, ¢, | dy fa) 
eens ede eee go ile i BOM ale! 


Wy 2182 0 YY 28 — 0.0888 
WY 21Y2 0 YiY1 21a — 0.1173 
WY 22 + 0.0132 YY eeu + 0.0366 
YY t12y + 0.0132 YY Ly Xq + 0.1040 
By 2 818» "+ 0.0166 ies oon + 0.1601 
2424 $1 Yo + 0.0443 fen bee 14a + 0.2135 
242, 8125 0 221 8124 + 0.1163 
212 Y1 85 — 0.0284 24s Y18q — 0.1423 
2% NY2 — 0.0494 4% Y1Va — 0.1087 
yey Y122 0 224 Y12%a — 0.1116 
2% 218 0 Bey Ziel — 0.0922 
Bye 214s 0 2424 21Ya — 0.1261 
22 212s + 0.0135 224 2124 + 0.0413 
22 UX» + 0.0128 224 %4Xq + 0.1012 
02, 818 + 0.0166 wyo, 818g + 0.1573 
U1, 81Y2 + 0.0443 112 81Ya + 0.2084 
Hy 8125 0 Dive 2324 + 0.1153 
U2, Y4 82 — 0.0284 0, Uy Yi 8a — 0.1532 
yey Y1Y2 — 0.0494 UX Y1Ya — 0.1054 
HX Y1 29 0 104 Y12%q — 0.1136 
2, 218 0 LyX, 218q — 0.0848 
LX, 21Ys 0 Hy Hy 21Ya — 0.1136 
112, 24 Za + 0.0128 D2 22a + 0.0371 
Uy Ht, X4Xy + 0.0135 Hy @42Xq + 0.1073 
ra US Se — 0.0002 K,8, Ky 8q — 0.0052 
3% Ki Ys — 0.0007 Ky 8, Ki Ya — 0.0101 
Ky, B25 0 K, 8; Ky 24 — 0.0056 
Bay; K, 8, — 0.0001 Kyy4 K, 8, — 0.0007 
Kyy, Ka Ys — 0.0002 Ky Yy1 Ki Ya — 0.0010 
aes Woaee 0 Ky Ei %, — 0.0008 
Ks ea Kes 0 Ka ea Ky 8, — 0.0004 
Siena Kiva 0 Ky% Ki Ya — 0.0008 
Ky Bey + 0.0000 Kies KAZ, + 0.0000 
Kya, Ki 2, + 0.0000 Ky, 2, Ky%q + 0.0005 
81 Y1 8182 — 0.0034 SY 8184 — 0.0293 
8141 81Y2 — 0.0078 8141 84a — 0.0262 
8141 812 0 8141 812a — 0.0248 
a1 Y18 + 0.0070 SY Y18a + 0.0484 
8141 Yi¥e2 +0.0151 841 Yi¥a ' + 0.0496 
8144 Y122 0 S1Y1 Y12a + 0.0432 
8144 28s 0 BY 21 8q + 0.0136 
811 21Y2 0 8141 21Ya + 0.0147 
8141 242 — 0.0019 © 84, 2124 — 0.0017 
811 X1X_ — 0.0019 81 Yy 12% — 0.0142 
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Table 23 (continued) 
ee SS SS eee eee 


Xb Xe Xa Xf (6, ¢1| ds fa) Xb Xe Xa re (6, 0; | dy fa) 


nnn EEE EEEEE EE SSEESESIES SS SSSSSSSSSSSS NEES 


8121 818, 0 8424 81 8q — 0.0162 
8421 81 Yo 0 812 81 Ya — 0.0248 
Bey 8125 + 0.0017 812 8124 + 0.0050 
82% Y1 82 0 812, Y1 8a + 0.0136 
$174 Y1Y2 0 $12 YiYa + 0.0147 
81,24 Y12 — 0.0019 812 Y1%a — 0.0017 
$124 2185 + 0.0020 8124 218q + 0.0314 
8124 Z1Y2 + 0.0070 8124 Ya + 0.0498 
8124 2125 0 8124 212 + 0.0118 
8121 Dies 0 8424 L,Xq — 0.0079 
Yy 24 818, 0 ies 81 8q + 0.0050 
911 51 Y2 0 Y1234 $1 Ya + 0.0062 
Yi2 8129 — 0.0006 Y124 82a — 0.0002 
Y1%1 Y182 0 Y1%1 18a — 0.0071 
i171 Yi¥2 0 Yi%1 Y1Ya — 0.0087 
Ye Yx 25 + 0.0010 Yy 24 Ys 2a + 0.0005 
Ys 24 2183 — 0.0006 Y424 21 8q, — 0.0078 
Yr 2 ZY ~— 0.0012 Yy 24 21Va — 0.0086 
Y12 2129 0 Yi 24 22a — 0.0030 
OPER Live 0 Ys ey 2 XL + 0.0022 
8,2 Sits + 0.0017 812, cafe, + 0.0187 
8,2, Ys — 0.0019 8,2, Y,Xq — 0.0142 
812, Pies 0 8,2, Bie — 0.0079 
812, L185 + 0.0020 812% X8q, + 0.0239 
812 L1Y2 + 0.0070 8,2, 11 Ya + 0.0373 
82 aoa 0 8,2, Dea + 0.0206 
YyXy 8X — 0.0006 YX, Sian — 0.0052 
YW, YX, + 0.0010 Yy XL, Ys La + 0.0078 
Yy Ly pa a 0 YX, Baie + 0.0023 
YX, L189 — 0.0006 YX, 18g — 0.0056 
YX, Ys — 0.0012 YyXy Ya — 0.0053 
YX, X22 0 YX, X12 — 0.0050 
2X, 81 Xp 0 2X, uaa — 0.0029 
@1Xy YX, 0 22, Yi %a + 0.0023 
24, 2a ite + 0.0004 Pane eae + 0.0039 
2, By X18, 0 Cai © 8g — 0.0031 
2X, X1Yo 0 2 Ly X4Ya — 0.0050 
Rie} Dies + 0.0004 2X4 Dita + 0.0009 
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